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Abstract 

We discuss the constraints on the time-varying equation of state for dark energy 
and the curvature of the universe using observations of type la supernovae from 
Riess et al. and the most recent Supernova Legacy Survey (SNLS), the baryon 
acoustic oscihation peak detected in the SDSS luminous red galaxy survey and cosmic 
microwave background. Due to the degeneracy among the parameters which describe 
the time dependence of the equation of state and the curvature of the universe, the 
constraints on them can be weakened when we try to constrain them simultaneously, 
in particular when we use a single observational data. However, we show that we 
can obtain relatively severe constraints when we use all data sets from observations 
above even if we consider the time-varying equation of state and do not assume a 
flat universe. We also found that the combined data set favors a flat universe even 
if we consider the time variation of dark energy equation of state. 



1 Introduction 



Almost all current cosmological observations indicate that the present universe is accel- 
erating. It can be explained by assuming that the universe is dominated by dark energy 
today. Although many candidates for dark energy have been proposed so far, we still do 
not know the nature yet. Many studies have been devoted to investigate dark energy as- 
suming or constructing a specific model and then study its consequences on cosmological 
observations such as cosmic microwave background (CMB), large scale structure (LSS), 
type la supernovae (SNela) and so on. On the other hand, many efforts have also been 
made to study dark energy in phenomenological way, i.e., as model independent as pos- 
sible. In such approaches, dark energy can be parameterized with its equation of state 
and constraints on it can be obtained using cosmological observations. Assuming that the 
equation of state for dark energy wx is constant in time, current observations give the 
constraint as wx ~ —1 [Ill21IHlllllSlin|- Although one of the most famous models for dark 
energy is the cosmological constant whose equation of state is wx = —1, most models 
proposed so far have time-varying equation of state. Thus, when we study dark energy, 
we should accommodate such time dependence in some way. 

Many recent works on dark energy investigate the time dependence of the dark energy 
equation of state using simple parameterizations such as wx = Wo + Wi{l — a) with a being 
the scale factor of the universe [ZIIHI. Parameterizing the dark energy equation of state in 
simple ways, the constraints on the time evolution of wx have been considered (for recent 
works on this issue, for example, see Ref. [Hj)- It should be mentioned that, when one 
studies the equation of state for dark energy, it is usually assumed that the universe is 
flat. It should be also noted that dark energy is usually assumed to be the cosmological 
constant when one derives the constraint on the curvature of the universe. However, it has 
been discussed that, assuming a non-fiat universe, the constraints on wx and the curvature 
of the universe can be relaxed to some extent even with the time independent wx from the 
CMB data alone ITU] . Also, even if we assume a fiat universe, there are degeneracies among 
the parameters which describe the evolution of dark energy, i.e., the time dependence of 
Wx when we consider the constraints on dark energy ^] . Furthermore it has been also 
discussed that if we remove the prior of a fiat universe, the degeneracies becomes much 
worse [ISj. Since it is very important to study the time dependence of wx to differentiate 
the models of dark energy and also the curvature of the universe to test the inflationary 
paradigm, we should investigate how the prior on the curvature of the universe affects the 
determination of the time-varying equation of state for dark energy and vice versa. Some 
works along this line have been done using a specific one-parameter parameterization for 
the time- varying wx [Ti] . 

In this paper, we study this issue, namely the determination of the evolution of dark 
energy and the curvature of the universe, in some detail using widely used parameterization 
of the equation of state for dark energy. We consider the constraints from observations 
of SNela reported in Refs. jH IH], the baryon acoustic oscillation detected in the SDSS 
luminous red galaxy survey [T3^ and recent CMB observations including WMAP . In the 
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next section, we summarize the analysis method we adopt in this paper. Then we discuss 
the constraints from above mentioned observations on the curvature of the universe with 
time-varying equation of state for dark energy followed by the analysis on wx without 
assuming a flat universe. The final section is devoted to the summary of the paper. 



2 Method 

In this section, we briefly summarize the method for constraining the parameters which 
describe the dark energy evolutions and other cosmological parameters. To study the evo- 
lution of dark energy, we use the following parameterization for the time- varying equation 
of state El 

Z 

wx{z) = Wo + ——wi = Wo + {1 - a)wi, (1) 

where z is the redshift. In this parameterization, the equation of state at the present time 
is Wxiz = 0) = Wq and for the early time it becomes wx{z = oo) = wq + Wi. Since we are 
interested in the late-time acceleration of the universe due to dark energy, we consider the 
case where the dark energy dominates the universe only at late time. Thus, in this paper, 
we assume 

Wo + Wi< 0, (2) 

in order not to include the possibilities of early-time dark energy domination. With this 
parameterization, the energy density of dark energy can be written as 



Pxiz) = pxoil + exp (^^) , (3) 



where pxo is the energy density of dark energy at present time. The Hubble parameter is 
given by 



H\z) = Hi 



+ zf + + zf + + zf + + exp [-^^^ 

where Ho is the Hubble parameter at the present epoch, f2j is the energy density of a 
component i normalized by the critical energy density and the subscripts r, m, k and X 
represent radiation, matter, the curvature of the universe and dark energy, respectively. 
To consider the constraints on dark energy and other cosmological parameters, we use the 
data from SNela, the baryon acoustic oscillation peak and the CMB. 

As for SNela data, we use the gold data set given in Ref. ^ and the first year data of 
the Supernova Legacy Survey (SNLS) released recently [Hj. Constraints from SNela can 
be obtained by fitting the distance modulus which is defined as 

M-m = Slog +25. (5) 
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Here c^l is the luminosity distance which is written as 



l + z 



(6) 



where S is defined as S{x) = sin(x) for a closed universe, S{x) = sinh(a;) for an open 
universe and S{x) = x with the factor a/P^ being removed for a fiat universe. 

We also use the baryon acoustic oscillation peak detected in the SDSS luminous red 
galaxy survey To obtain the constraint, we make use of the parameter A which is 
defined as 



A 



-.s 



dz' 



H{z')/H, 



2/3 



(7) 



where zi = 0.35 and A is measured to be A = 0.469 ± 0.017 [TK] . 

For the CMB data, we only use the shift parameter R which determines the whole shift 
of the CMB angular power spectrum fTBj. R is given by 



R 



where Z2 = 1089. It has been discussed that using the CMB shift parameter is a robust way 
to include the constraints from observations of CMB ^^I- From the recent observations 
of CMB including WMAP, CBI and ACBAR, the shift parameter is constrained to be 
R= 1.716 ±0.062 Pini 

In this paper, we only consider the effect of the modification of the background evolu- 
tion by the change in the cosmo logical parameters. In fact, the properties of dark energy 
can also modify the evolutions of cosmic density fluctuation. When we consider the effects 
of dark energy perturbation, we also have to specify the speed of sound of dark energy 
component. However, such modiflcation can arise at low multipole region of the CMB 
power spectrum where the errors due to cosmic variance are large. Thus the constraints 
from CMB on dark energy mostly come from the position of acoustic peaks which can 
be described by the shift parameter. Hence we do not consider the perturbation of dark 
energy in this paper*^ . 



3 Constraints on the vs. Qx plane 

Now we discuss the implication of the dark energy evolution on the determination of the 
curvature of the universe. For this purpose, we derive the constraints from the observations 
on the Qrn vs. Qx plane. First we consider the case with the constant equation of state for 

T^^However, the speed of sound can be useful to differentiate the models of dark energy. There are some 
works which discuss the constraint on the speed of sound. For interested reader, we refer Refs. [ TSIIT^I^ . 
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Figure 1: Constraints on Qrn and fix from SNela (a), the baryon acoustic oscillation peak 
(b), the CMB shift parameter (c) and all data combined (d). For the constraint from 
SNela data and all data combined, we show the contours obtained from the gold set of 
Ref. Ill and SNLS |n| separately. Contours of la (red solid line) and 2a (green dashed line) 
are shown (for SNLS data, la and 2a contours are shown in blue dash-dotted line and 
purple dotted line, respectively). We assumed the cosmological constant as dark energy 
and a flat universe here. 



dark energy. In Fig. ^ we show the contours of la and 2a constraints from observations of 
SNela (a), the baryon acoustic oscillation peak (b) and the CMB shift parameter (c). We 
also show the constraint from the combination of all data sets (d) . In Fig. Q we assume 
the equation of state as wx = — 1- For the constraints from SNela, we used the gold data 
set from Ref. |4j and the data from SNLS 6j separately. As we can see, although each 
constraint from a single observational data has the degeneracy in the Qm vs. Qx plane, 
when all the data sets are combined, we can obtain a severe constraint as Qm ~ 0.3 with 
fifc ~ 0, which is a well-known result. 

Next we consider the case where the equation of state for dark energy is allowed to 
vary, but still we keep wx constant in time. In Fig. |21 we show the contours of la and 2a 
constraints after marginalizing over the values of wx- Here we assumed the prior on wx 
as — 5 < Wx < 0. As seen from the figure, the allowed regions become larger compared 
to those for the case with wx = — 1 when we use a single data set alone. However, by 
combining all data, we can obtain almost the same constraint as the case with wx = —1 
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Figure 2: The same as Fig. ^except that we marginahzed over the value of wx- Here we 
considered the constant wx- 

even though we marginahze over wx- This can be understood by noting that the values 
of Wx which give the minimum for fixed (fim, ^x) from each data set are different. To 
see this, we also plot contours of constant wx which gives minimum values of at each 
point on the flm vs. Qx plane in Fig. El We can clearly see that the favored values of wx 
vary for different data sets. Observations we consider here measure some distance scales 
to certain redshifts which are determined by the energy density of matter Q^n, dark energy 
Qx and the equation of state wx- If we vary the value of wx, the density parameters flm 
and Qx can have more freedom to be consistent with observations since the fit to the data 
depends on the combinations of these quantities. For a single observation, when wx is 
allowed to vary, larger range of {flm, f^x) can be consistent with observations. However, 
when we use all observations, the combinations of flm, f^x and wx which are consistent 
with observations become fairly limited. Thus while the allowed regions become larger 
for each data set, when we combine all the data, the allowed region converges towards 
the concordance model with flm + fix ~ 1 and flm ~ 0.3 even if we do not assume the 
cosmological constant as dark energy. 

Next we discuss the case with the time dependent equation of state parameterized as 
Eq. (d}. As it has been already pointed out in the literature, even if we assume a flat 
universe, there exists a degeneracy among the parameters which describe the evolution of 
equation of state for dark energy using a single observational data. It is also known that 
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Figure 3: Contours of constant Wx which gives minimum when we marginahze over 
Wx for the constraints from SNela (a), the baryon acoustic oscillation peak (b), the CMB 
shift parameter (c) and all data combined (d). Contours of wx = —2 (red solid line), — 1 
(green dashed line) and —0.8 (blue dotted line) are shown except for the case where all 
three data are combined (panel (d)). For the panel (d), contours of wx = — 1.1 (red solid 
line), —1 (green dashed line) and —0.9 (blue dotted line) are shown. For SNela, the data 
from SNLS are used here. 



the degeneracy can be removed using more than one observation assuming the flatness. 
Here we discuss to what extent the evolution of dark energy equation of state can affect 
the determination of Qm and Qx when we do not assume a flat universe. In Fig. |3J the 
contours of la and 2a constraints are shown as the same manner as Fig. ^ except that 
we varied the both values of Wq and Wi which appear in Eq. (P) and marginalized over 
them to obtain the constraint. We assumed the prior on them as —5 < < and 
—4 < Wi < 4 under the condition of Eq. For the constraint from SNela, the allowed 
regions get larger compared to the case with the constant equation of state discussed 
above (see Fig. |21). As for the constraint from the baryon acoustic oscillation peak and the 
CMB shift parameter, the allowed regions are almost the same as the case with a constant 
equation of state. For observations which measure a single distance scale, the fit to the 
data does not significantly become better, even if the time evolution of wx is allowed, 
due to the degeneracy between wq and wi. Furthermore, when we use all three different 
observations together, the allowed region becomes almost the same as that of the case 
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Figure 4: The same as Fig.Hlexcept that we consider the time dependent equation of state 
for dark energy and marginahze over Wq and Wi. 



with a constant equation of state. This is because each observation favors different values 
of wx- Even if we consider the time- varying equation of state, different combinations of 
Wo and wi are chosen to minimize the value of for each observation. Thus the allowed 
region from the combined data set is almost unchanged although the constraints from a 
single observation, in particular from SNela, can become weaker. 

Here we comment on the effect of the prior on wq and Wi. Although the constraints 
from a single data alone somewhat depend on the prior, the allowed region for all data 
combined is almost unaffected since the combinations of wo and wi favored around the 
allowed region from all data sets are far from the the edge of the prior we assumed. 



4 Constraints on the vs. wx plane 

In this section, we consider the constraints on the Q^n vs. wx plane. First we show the 
constraint on Qrn and Wq in a flat universe with the equation of state for dark energy being 
constant. In Fig. El contours of la and 2a constraints from observations of SNela (a), the 
baryon acoustic oscillation peak (b), the CMB shift parameter (c) and all data combined 
(d) are shown. As is well-known, SNela and CMB are complementary for constraining 
dark energy, which can be seen from the figure. In addition, we can see that the constraint 
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Figure 5: Constraints on and wx from observations of SNela (a), the baryon acoustic 
oscillation peak (b), CMB (c) and all data combined (d). We assumed a flat universe and 
a constant equation of state here. 

from the baryon acoustic oscillation peak is also complementary. Thus we can obtain a 
severe constraint using all three data sets. 

Next we consider the case with the time-varying equation of state for dark energy. 
Here we discuss the case with a fiat universe. In Fig. IHl the constraints on and Wq 
are shown marginalizing over Wi. When we marginalize the values of Wi, we assumed the 
prior on Wi as —4 < Wi < 4. Similarly to the situation where we constrain fl^n and Qx 
discussed in the previous section, if we consider the time- varying equation of state for dark 
energy, the allowed region becomes larger when we use a single data set. However, when 
we consider all three data sets, we can obtain a relatively severe constraint. Although 
the allowed region of wq extends to larger values, flm is constrained to be flm ~ 0.3 
which is the almost same as the case with assuming a constant equation of state. This 
is because the values of wi favored by different data sets given and wq are different. 
Thus the degeneracy among the parameters describing the time dependence of the dark 
energy equation of state can be removed drastically by using all three data sets. 

Now we discuss the constraints without assuming a fiat universe. First we discuss the 
case with wx being constant in time. In Fig. [71 we show the constraints from SNela (a), 
the baryon acoustic oscillation (b) and all data combined (c). For the last case, the data 
from the CMB shift parameter is included. To obtain the constraint, we marginalized 
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Figure 6: The same as Fig.|31except that we consider the time dependent equation of state 
for dark energy as in Eq. (^. To obtain the constraint, we marginahzed over Wi. Here we 
assumed a flat universe. 

over Qk with the prior —0.3 < fi^ < 0.3. Here we do not show the constraint from the 
CMB shift parameter alone since we cannot obtain a significant constraint from it in the 
parameter region we consider in Fig. [7| This is because, as it has already been pointed 
out, the curvature of the universe and wx are strongly degenerate in the CMB power 
spectrum jTU]. Again, although the constraints from a single data set alone are weakened, 
if we consider all data sets, we can obtain almost the same constraint as that with the 
case where a fiat universe is assumed. Notice that the constraint on wq is also not changed 
much compared to that with flat universe prior. Hence from Fig. we can say that the 
prior on the curvature does not affect the determination of the equation of state for dark 
energy much. 

We also show the case where we consider the time-varying equation of state without 
assuming a flat universe. In Fig. |H1 the constraints are shown as the same as Fig. [7| except 
that we marginalized over and Wi in this case. The prior on wi is taken as —4 < Wi < A. 
The CMB shift parameter cannot give a signiflcant constraint in the parameter range we 
consider in this case too. As in the previous cases, using each observational data alone, 
the allowed region becomes signiflcantly larger compared to those with constant equation 
of state and a flat universe. However, when we use all data sets, the allowed region does 
not change much. From Figs. El and we can conclude that the prior on the curvature 
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Figure 7: Constraints on Q^n and Wx from observations of SNela (a), the baryon oscillation 
(b) and all data combined (c). Here we do not assume a flat universe, but a constant 
equation of state for dark energy is assumed. In this figure, we marginalized over Qj^. 



of the universe does not affect the constraint on fi^ and wq much if we use all data sets. 
In particular, we can obtain the constraint such that Qm ~ 0.3 without assuming a fiat 
universe. Moreover, the assumption of the constancy of the equation of state for dark 
energy also does not affect the constraint on Qm- On the constraint on the equation of 
state, if we consider the time- varying equation of state, the allowed region becomes slightly 
larger even if all data combined is used. 



5 Constraints on the evolution of dark energy 

In this section, we consider the constraints on the time dependence of dark energy equation 
of state, i.e., on wo and wi without assuming the prior of a fiat universe. However first 
we discuss the case with a fiat universe for the later comparison with the case where a 
flat universe is not assumed. In Fig. El contours of la and 2a allowed regions are shown 
for the case with a flat universe. Here we fix the energy density of matter as Qm = 0.28. 
Notice that we do not consider the region where Wq + Wi > in order not to include 
the possibilities of early-time dark energy domination. As we can see from the figure, the 
constraint from SNela is stringent compared to baryon acoustic oscillation peak and CMB. 
Since the position of the baryon acoustic oscillation peak measures a single distance scale 
from z = to z = 0.35 and the CMB shift parameter also gives a single scale from 2; = to 
z = 1089, there is strong degeneracy among the parameters which describe the equation 
of state for dark energy. Notice that the distance scales are determined by the integration 
of the inverse of the Hubble parameter, which can smear out the information on the time 
dependence of the equation of state. This is the reason why the strong degeneracy exists 
when only a single scale is considered. However, as for SNela data, the degeneracy is 
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Figure 8: The same as Fig.[7|except that we consider the time dependent equation of state 
for dark energy and do not assume a flat universe. In this figure, we marginahzed over Qj^ 
and Wi. 



removed to some extent since we have the distances from z = to various redshifts*^ . 
Thus observations of SNela can mainly constrain the time dependence of the equation of 
state for dark energy. 

Next we show the constraints on wq and Wi in a flat universe without assuming a 
particular value for ^2^- To obtain the constraint, we marginalized over with the prior 
< Qrn < 0.5. In Fig. EH we show the constraints on wq and Wi from SNela (a) and all 
data combined (b). We do not show the constraints from the baryon acoustic oscillation 
peak and the CMB shift parameter because we cannot obtain significant constraints from 
them in the parameter region of Wq and Wi we consider here. Also in this case, such severe 
degeneracies among the parameters exist because those observations measure a single 
distance scale. As in the previous case, when we include all data sets in the analysis, 
we can obtain almost the same constraint as that with being fixed to 0.28. This is 
partly because the best fit value of flm is near 0.28 when we marginalized over it, but 
also because the values of Qm minimizing for each observational data set are different. 
To see this, in Fig. ^2 we show contours of constant Qm which minimize when we 
marginalize over for the constraints from SNela (a), the baryon acoustic oscillation 
peak (b), the CMB shift parameter (c) and all data combined (d). As is clear from the 
figure, preferred values of flrn for each observation are quite different. Thus when all data 
are used, we can obtain a severe constraint although a single observation cannot constrain 
the equation of state for dark energy. It should be also mentioned that the prior on Q^n 
does not affect the constraint on the time dependence of the equation of state when we 
use all data combined. 

Now we discuss the case with Qrn not being fixed and not assuming a flat universe. In 
Fig. ^[ the allowed regions are shown after marginalizing over Qm and Qk- We assumed 

^^Of course, there is a limitation to determine the time dependence of the equation of state using SNela 
data lllllig. 
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Figure 9: Constraints on wq and Wi from SNela (a), the baryon acoustic oscillation peak 
(b), the CMB shift parameter (c) and all data combined (d). We assumed flm = 0.28 with 
a fiat universe. 

the prior on these variables as < Qm < 0.5 and —0.3 < < 0.3. In this case, the 
constraint is significantly weakened when we consider a single data set alone. We only 
report here the constraint from SNela (a) and that from all data combined (b) since the 
baryon acoustic oscillation peak and the CMB shift parameter cannot give meaningful 
constraints on this plane in this case too. As we can see, when we consider the constraint 
from all data sets, it significantly becomes severe compared to that from SNela alone as 
seen from the figure. This is because the favored values of Qm and from each data 
set are different as in the previous cases. Thus we can conclude that the prior on the 
curvature does not affect much the determination of the equation of state for dark energy 
if we use all data combined. 

6 Summary 

We considered the constraint on the curvature of the universe and the equation of state for 
dark energy from observations of SNela, the baryon acoustic oscillation peak and the CMB 
shift parameter. Usually, when one discusses the curvature of the universe, dark energy is 
assumed to be the cosmological constant. Moreover, when one considers the constraint on 
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Figure 10: Constraints on wq and Wi from observations of SNela (a) and all data combined 
(b). A flat universe is assumed but we marginalized over Q^- 



the evolution of dark energy, in particular the time dependence of the dark energy equation 
of state, a flat universe is usually assumed. In this paper, we discussed the constraints 
on the curvature of the universe without assuming the cosmological constant and also the 
time dependence of the equation of state for dark energy without assuming a flat universe. 
We showed that the constraint on the curvature of the universe is significantly relaxed 
from a single observation when we allow the time dependence in the dark energy equation 
of state. However, it was also shown that, when we use all data sets, the curvature of the 
universe or the energy density of matter and dark energy are severely constrained to be 
Qrn ~ 0.3 with a flat universe even if we consider the time- varying equation of state for dark 
energy. Observations we consider here measure some distance scales to certain redshifts 
which can be determined by the energy density of matter Qrn, dark energy Qx and the 
equation of state wx- If we assume a broad range for wx, the energy density of matter Qrn 
and dark energy Qx can have more freedom to be consistent with observations since the 
fit to the data depends on the combinations of these variables. For a single observation, 
if Wx is allowed to vary more freely, much larger range of values of Qm and Qx can be 
consistent with observations. However, when we use all observations, the combinations 
of Qm, and wx consistent with observations become fairly limited. Thus even if the 
allowed regions become large for each data set, when we combine all the data, we can 
obtain a severe constraint, which is interestingly almost the same region as that we obtain 
assuming the cosmological constant as dark energy. 

We also investigated the constraint on the time-varying equation of state for dark 
energy without assuming a fiat universe. Similarly to the situation where we constrained 
the curvature with the time-varying equation of state, the allowed region for wx becomes 
larger when we use a single observational data. However, if we use all data sets considered 
in this paper, we can obtain almost the same constraint as that in the case where a fiat 
universe is assumed. 

Finally, we summarize what we found in this paper. The combination of the current 
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Figure 11: Contours of constant Qm which gives minimum when we marginahze over 
Qm for the constraints from SNela (a), the baryon acoustic oscillation peak (b), CMB shift 
parameter (c) and all data combined (d). Contours of flm = 0.3 (red solid line), 0.2 (green 
dashed line) and 0.1 (blue dotted line) are shown except for the case where all three data 
are combined (panel (d)). For the panel (d), flm = 0.29 (red solid line), 0.28 (green dashed 
line) and 0.27 (blue dotted line) are shown. For SNela, we use the data from SNLS. 

observations 

• favors a flat universe regardless of the prior on the equation of state for dark energy 
with or without time evolution. 

• favors Qm ~ 0.3 regardless of the flatness prior and the prior on the equation of state 
for dark energy with or without time evolution. 

• yields constraints on the time evolution of dark energy equation of state regardless 
of the prior on Qrn and Qk- 

In future observations, we can obtain much more stringent constraints on the equation 
of state for dark energy as well as the curvature of the universe regardless of the prior on 
other cosmological parameters. Hence we can expect that we will be able to have much 
insight on the nature of dark energy and also the inflationary paradigm. 
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Figure 12: The same as Fig. except that we marginahzed over and Q^- 
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